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Abstract. Lie group analysis of differential equations is a generally 
recognized method, which provides invariant solutions, integrability, 
conservation laws etc. In this paper we present three characteristic 
examples of the construction of invariant difference equations and 
meshes, where the original continuous symmetries are preserved in 
discrete models. Conservation of symmetries in difference modeling 
helps to retain qualitative properties of the differential equations in 
their difference counterparts. 

1. Introduction 

Symmetries are intrinsic and fundamental features of the differential equations 
of mathematical physics. Consequently, they should be retained when discrete 
analogs of such equations are constructed. 

The group properties of a heat transfer equation with a source 

u t = (K(u)u x ) x + Q(u), (1) 

were considered in [4], and all choices of K(u) and Q(u) which extend the sym- 
metry group admitted by the general case of equation (1) were identified. In this 
paper we consider two partial cases of nonlinearities from the classification in [4] 

u t = (u a u x ) x ± u n , a, n = const, (2) 

u t = u xx + Suliau, 5 = ±1, (3) 

together with linear case 

u t = u xx , (4) 
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whose group properties were known by S. Lie. For all cases we construct dif- 
ference equations and meshes (lattices) that admit the same Lie group of point 
transformations as their continuous limits. 

We recall that Lie point symmetries yield a number of useful properties of 
differential equations [13,16,10]: 

a) A group action transforms the complete set of solutions into itself; so it is 
possible to obtain new solutions from a given one. 

b) There exists a standard procedure to obtain the whole set of invariants and 
differential invariants for a symmetry group; it yields the invariant representation 
of the differential equations and the forms of invariant solutions in which they 
could be found (symmetry reduction of PDE). 

c) For ODEs the known symmetry yields the reduction of the order; if the 
dimension of symmetry is equal to (or greater then) the order of ODE, then we 
have a complete integrability. 

d) The invariance of PDEs is a necessary condition for the application of 
Noether's theorem on variational problems to obtain conservation laws. 

e) It should be noticed that Lie point transformations have a clear geometrical 
interpretation and one can construct the orbits of a group in a finite dimensional 
space of independent and dependent variables. 

The structure of the admitted group essentially effects the construction of 
equations and grids. Group transformations can break the geometric structure of 
the mesh that influences the approximation and other properties of a difference 
equation. Early contributions to the construction of the difference grids based on 
the symmetries of the initial difference model are [6,8]. Classes of transformations 
that conserve uniformity, orthogonality, and other properties of meshes will be 
defined below. 

In accordance with equation (JIJ) we consider Lie point transformations in a 
space with two independent variables: t and x. Let 

be an operator of a one-parameter transformation group. Dots denote prolon- 
gation of the operator on other variables used in the given differential equation: 

F(x, t, u t , u, u x , u xx ) = 0. (6) 

The group generated by (JHJ transforms a point (x,t,u,u t ,u x ,u xx ) to a new one 
(x*, t*, u*, ul, u* x ,u xx ) together with equation (0). This situation changes when 
applying Lie points transformations to difference equations. Let 

F(z) = (7) 

be a difference equation defined on some finite set of points z 1 , z 2 , . . . (difference 
stencil) on a mesh. In contrast to the point (x,t,u,u t ,u x ,u xx ), the "difference 
point" - the difference stencil has its own geometrical structure. 
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Let 

tt(z,h)=0 (8) 

be an equation that define a difference stencil and a mesh. As an example it will 
be a uniform mesh if the left step (spacing) equals the right step: 

h + = hr. (9) 

The invariance of the difference equation (jJJ) depends on the invariance of (jHJ), 
since the latter must be included in the general condition of invariance: 

XF ^\ 7m = 0, 
XCl{z,h) ]mw = 0. [W) 

Relations between the two conditions in (fTUJ) depends on whether Q, £*, £ x 
depend on solution or not. If Q u = = 0, then the conditions (fTUJ) could 

be considered independently. 

Thus, what makes our approach [5-8] special is the inclusion of the second 
equation of (jlOj) in the conditions of invariance, admitting the whole set of prop- 
erties a)-e) stated for equations ((Zj) and (JHJ). 

There exist a few ways to avoid transformations of the difference stencil and, 
consequently, transformations of a mesh. One way is to restrict transformations 
to the case when independent variables are not changed: = £ x — 0, yielding any 
mesh invariant (see [12]). But this restriction is very strong and would exclude 
most symmetries of physical problems. 

Another approach is connected to evolutionary vector fields. It is known 
[10,16] that the symmetry operator could be represented as 

x={n-eu t -? s u x )^ + '-- (ii) 

which could be viewed as representative of a factor algebra by the ideal 

e(z)D t + C(z)D x , (12) 

d d d d 

where D t = — + Ut— — V ■ ■ ■ , D x = — — h u x — . . ., both admitted by all equa- 

ot ou ox ou 

tions (JBJ). For differential equations this approach gives an equivalent result but 

with some losses in points a)-e). In particular we lose the geometric sense even 
for point transformations due to the realization of group transformations in in- 
finite dimensional spaces (t,x,u t ,u x ,u tx ,u xx ,u tt , . . .), there is no procedure for 
calculating invariants of operators of the type given by (jTTjl . etc. 

As it was shown in [5], the exact representation of the operator of total differ- 
entiation D x in a space of difference variables is given by the following operators 
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(the same is also true for D t ): 



(13) 



ox +h Ou +h / — ' n +h 

n=l 

£- = — + £(«)— + •••, £ = V D n , 

Ox -h ou -h ^— ' n -h 

n=l 

where D and £) are right and left difference derivative operators on a uniform 

+h -h 

mesh. As difference operators are related to corresponding shift operators S , S '■ 

+h -h 

(s-i) 

±h n 

we obtain another representation for D: 

±h 

d = E ^ n = E - z} i—(s -ir, (14) 

n=l n=l 

oo (1- S)™ 

n=l 



The group transformation for operators (|13jl can be obtained by the exponen- 
tial mapping or by means of the so called Newton series (see [5]). It is important 
to notice that every difference equation on a regular mesh admits the operators 
D + and D~, which do not change a mesh. It was shown [5], that a family of 
operators £(z)D ± forms an ideal in the Lie algebra of operators (j3J), since it is 
possible to rewrite these operators as evolutionary vector fields 

X=(v~ ^Di{u) - CDt(u)) !; + •••■ (16) 

(We follow here the right semiaxis representation D + . The same formulation can 
be done with help of left semiaxis representation D~ .) 

It is important to notice that the representation (|13|) is true only for regular 
(uniform) meshes. Thus, the evolutionary vector fields are applicable only 
for the groups which do not change uniformity of a mesh. So, one can not apply 
them for modern numerical methods with moving meshes, self adaptive meshes 
or multi-grid methods etc. 

Let us consider a one-parameter transformation group which is generated by 
the operator 

x = ^l + -- < 17 > 
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We prolong (JT7j) on the right and left steps h + , h by means of relations h + = 
x + — x and h~ = x — x~ , where f^=S (/): 

±h 

x -4 x + 4u + - + &«> - *m + 1? " 5« )] 5F' (18) 

From (J 18)) it is easy to obtain the invariance condition for a uniform mesh in 
a given direction. Let be invariant with respect to (jl8|) . then 

5(0" 2^+5(0 = or £>£>(£) = 0. (19) 

+h -h +h-h 

Condition (|19p is a strong limitation on the admitted group. In addition the 
coefficients of (J16j) are the power series of D or S , since one should consider the 

±h ±h 

whole set of mesh points and not a stencil only. 

Let us illustrate the above with a simple example. The ODE 

u xx = it 2 , (20) 

can be viewed as the stationary case of equation (0) with a = 0, n = 2. Equation 
(120)) has the following Lie point symmetries: 

Xi = P, X 2 = x^--2u^-. (21) 
ox ox Ou 



As a difference analog of (J20)) we consider 

u 2 . (22) 



u + — 2u + u 2 



h 2 

on a uniform mesh 



(23) 



where u + = S(u), u = S{u). 

+h -h 

The equations ()2^j) and (j2Hl) use a three-point stencil or subspace (x, x + ,x~, u, u + , u~ 
and the operators (|21)) have the following prolongation for the shifted points of 
the difference stencil: 



d d d d , d d 

X\ = 7T- + TT-T + ~ i X 2 = x— + x^—— + x 



dx dx + dx~ ~ dx dx + dx~ 

-2u^--2u+J--2u-J- + h+^ + h- 9 



(24) 



du du + du~ dh + dhr 



The symmetry algebra (|2Tj) acts on the space (t, x), so the coefficients of (|24j) 
have the same form in different points of the stencil. The prolongation forms for 
h + and h~ are easily derived from the relations h + = x + — x, h~ = x~ — x. 
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It is easy to verify that the equations (|22|) and (|23|) are invariant under the 
action of 

— 2u + M~ , 
A- | - u 2 \ = 0. 



X 2 {h+-h-) ]trm = 



(25) 

= n 

1221 



(Operator X x leaves (j2^|) and unchanged). 

It follows that the system (J22|) and (|23|) has the same Lie point symmetry 
as its continuous limit. Notice that the invariance conditions ()25|) are mutually 
independent. 

As for the continuous case one can easily calculate the finite difference invari- 
ants for (J24j) by solving the standard system: 

Xi I{x,h + ,h-,u,u + ,u-) = 0, i = 1,2. (26) 

The solution of (|26|) yields the whole set of difference invariants 

1' = ^, I 2 = ~, I 3 = ~, I* = (h+)\ (27) 
a u u 



It follows that the difference model ()22|) and ()23|) can be represented by means 
of the invariants (JZZD as I 2 + I 3 - 2 = I 4 and J 1 = 1. 

Let us now consider the evolutionary vector fields for the difference equation 
(I22|) (we consider the right semiaxis representation D + ). We prolong the operator 
(fT3"j) on all points of a given stencil (x, x + , x~, u, u + , u~): 

D + — ^ + ^ + ^ (28) 

<9x dx + dx~ x du x 0u + x 0u~ ' 

where 

n=l n=l n=l 

The evolution vector fields (J 16)) will have the following forms: 

- 

X x = -Xx + D + = u x — + u+ — — + w~ 



h" 11 x du x du+ x du-' 

X 2 = -X 2 + xD+ = -h + ^--h-^- (29) 
h 0h + 0h~ 



+ (2u + xu x )— + (2u + + xu x )^— + (2u~ + xu~ s 



Ou y x '0u+ v x '0u 
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It is not easy to check the invariance conditions of equation for the 
operators because one should use not only the equation but all its 

sequences obtained by shifting to the right. A harder question is how to produce 
the finite difference invariants (J27J) by means of (|29|1. That is why we prefer the 
first classical way for Lie point symmetries and include a mesh in the invariance 
condition (JTUJl (a very similar approach for the semi-discretized nonlinear heat 
equation was introduced recently in [3]). An additional convincing reason to apply 
the classical representation of Lie point symmetries springs from comparison of 
unsuitable evolutionary vector fields approaches to invariant variational problems, 
developed in [5], and a clear classical way to construct Noether type theorems 
for difference equations [7]. 

Another approach to the symmetry of difference equations on a fixed uni- 
form mesh was introduced in [9]. However, that way is only applicable to linear 
equations; moreover, it requires knowing a complete set of solutions of difference 
equations. The newly introduced approach [11] deals with evolutionary vector 
fields on a uniform mesh. The advantage of the last two approaches seems to be 
in the potential of finding non-point symmetries of difference equations which are 
not available in their continuous limits. 

Returning to equation ((H), we notice that a transformation defined by (JHJ) 
conserves uniformity of a grid in t and x directions, if 

= 0, (30) 

+T — T 



DD{e) = 0, (31) 

+h-h 

where D and D are right and left difference operators in t direction. 

Conditions ()30|) and ()31j) are not sufficient to describe the invariance of an 
orthogonal mesh. For an orthogonal mesh to conserve its orthogonality under 
the transformation, it is necessary and sufficient that [6,8]: 

D(e) = -D(0- (32) 

+h +t 

When condition (J32j) is not satisfied for a given group, the flatness of the 
layer of a grid in some direction is rather important. For evolution equations it is 
significant to have flat time layers, since otherwise, after transformations, some 
domains of a space could be in the future, while others in the past. We have a 
simple criterion [6,8] for the preserving the flatness of the layer of a grid in the 
time direction under the action of a given operator ©: 

= 0. (33) 

+h+T 

So, the conditions (J5U|l - (J53j) provide a geometry of grids that is based on the 
Lie group symmetry. These conditions will be used in what follows. 
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2. Invariant model for the equation Ut = {u a u x ) x ± u 11 

The equation 

u t = (u a u x ) x ± u n , <j,n = const, (34) 

admits a 3-parameter symmetry group. This group can be represented by the 
following infinitesimal operators [4]: 

x '=l ; x >=h *= 2 ("- i '4 + ("- ff - i) 4- 2 4- < 35 > 

The set (|33j) satisfies all conditions (j3llj) - (j3^j) . So, we can use an orthogonal grid 
that is uniform in the t and x directions. Let us consider the set of operators (|35p 
in the space (t,t,x, h + , h~ ,u,u + ,u~,u,u + ,U-) that corresponds to the stencil 
shown in Fig. 1. 

(x — h,i,u_) (x,i,u) (x + h,i,u + ) 



(x — h,t,U-) (x,t,u) (x + h,t,u + ) 

Figure 1: The stencil of the orthogonal mesh. 
There are 7 difference invariants of the Lie algebra (|35l) : 



I^; ™~; «; % ^; St; U -. (36) 
n u u u u u 

The small number of symmetry operators (|35|) provides us with a large number 
of difference invariants (|36p. Thus we are left with some additional degrees of 
freedom in invariant difference modeling of ()34j) . By means of the invariants ()36|) . 
we could write the following explicit scheme for (J33j) : 

' ' 1 Mx - — =— «s ±« , (37) 



where u 



r h\\ 2 J h \ 2 J h 

u+ — U U — U- 



h h h h 

This scheme is not unique and one could construct another form of invariant 
difference equation. For example an implicit scheme could be as follows: 

^—1 = J_(^+i _ 2^+1 + u°_ +1 ) + u n (38) 
t n z 
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Notice, that continuous limit of the last difference equation 

u t = {u° +1 ) xx + u n . (39) 

is equivalent to the equation (J%4"|) up to the scaling of x. But scheme (13%)) is 
not equivalent to the scheme ()37|h because there is no point transformation that 
relates them. In [18] Samarskii et all considered the case n = a + 1, a > and 
found a finite-difference blow-up solution for the equation ((3%J that is invariant 
with respect to the operator 

where T is constant. The operator (J4Uj) defines a subgroup which is equivalent 
to the self-similar subgroup with the operator X 3 of the set (J33J). Let us find the 
solution of the problem ()38|) in the invariant form: 

1 

»=(l"|) 

This solution is sought ([18]) on the time mesh with the infinite number of de- 
creasing steps 

T * = a ^ p ~ aj > J = °>M.-> ( 41 ) 

where p > 1 is constant. This mesh corresponds to the blow-up time 

T = ° p ~ 1 
p p CT - 1 ' 

For the function 6(x) we have the equation 

(Oxx + ^ +1 = -^- ( 42 ) 



The solution of the problem ()42|) was found in [18] for the case a = 2. Let us fix 
an arbitrary I 
length equals 



an arbitrary M > and let h — 2 sin I . ) . In this case the localization 

\2(M+1)/ 



, 37rh ( ■ n 

i/, = arcsm — , U < h < 2 

2 V 2/ 



(see [18]). Then, one can verify that the solution of the problem ()42|) in the points 
Xk = kh has the form 

fc/l = V2^1-^sin 2 ^JJ sm(a h kh), k = 0, 1, M + 1, (43) 



where ah = ir/lh- 

The obtained function u gives the blow-up solution of the problem in the case 
a = 2, n = 3, I = lh- This solution tends to infinity in all points of the space 
grid, conserving the structure. As h — > 0, the difference solution ()43)) tends to 
the solution of the ordinary differential equation: 

/g\ 1/2 

e(x)=l-j sin (|), < a; < / = 3tt. 



3. Invariant difference model for the semilinear heat transfer equation 

The semilinear heat transfer equation 

u t = u xx + 5u\iau, 6 — ±1, (44) 

admits the 4-parameter Lie symmetry group of point transformations [4], that 
corresponds to the following infinitesimal operators: 

x « ; * = 2e ft» ,«« » i, = Ai (45) 

at ax ax aw aw 

Before constructing a difference equation and a grid that approximate (pH^ and 
inherit the whole Lie algebra ([45)1 . we should first check condition (}3*2*|) for the 
invariance of orthogonality. The operators Xi, X2 and X4 conserve the orthogo- 
nality, but X 3 does not: condition is not true for the last operator. Conse- 
quently an orthogonal mesh cannot be used for the invariant modeling of (J44I) . 

The conditions (jHB^l are true for the complete set ()45|). so it is possible to use 
a nonorthogonal grid with flat time layers. An example of a grid with with flat 
time layer is shown in Fig. 2. 




x 



Figure 2: An evolutionary mesh with flat time-layers. 
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A possible reformulation of equation (|44|) by using the four differential invari- 
ants in the subspace (t dt, dx, du): 



u J u u dt udt 4 V dt I 



is given by the system: 



or 



dx 
~dt 

du 
~dt 



J 3 
J 1 



0; 



-2- 



It 



u 

u xx + 5u\n.u — 2—. 

u 



(46) 



So, the structure of the group (|43j) suggests the use of two evolution equations. 

As the next step, we will find difference invariants for the set X1—X4 of the 
group (J45)) . These invariants are necessary for the approximation of the system 
). We will use the 6-point difference stencil, as shown in Fig. 3 



(x — h ,t,u_) (x,i,u) (x + h + ,i,u + ) 



(x — h ,t, U-) (x, t, u) 



[x + 



U+j 



Figure 3: The stencil of the evolutionary mesh. 

on which we will approximate the system (J46j) . The stencil defines the difference 
subspace (t, t, x, x, h + , h~ , h + , h~ , u,u + ,u_, u,u + ,uJ) . The group (|45|) has the fol- 
lowing difference invariants in this subspace: 

I 1 = t; I 2 = h + ; I 3 = h~; J 4 = h + ; J 5 = h~; 

I 6 = (lnu) x - (ln«) s ; I 7 = {lnu) x - (ln«) s ; 
t8_xa. , o/ir h~ ... , h+ 



/- = ^ + 2(e--l)^ FTF (]n«). + FTF (]n«) I 

I 9 = 5Ax + 2(1 - e~ 5T ) I , - s (\nu) x + s (Inn) 

+ h+ + h- 

I 10 = 5(Ax) 2 + 4(1 - e~ 5r ) (in u - e 5r In w) ; 
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\nu + — \nu Inu — lnu- 
where Ax = x — x, (mu) x = — , (mu)x = — . 

To obtain an invariant difference model, it is natural to use the difference invari- 
ants. An explicit model is given by 

I 8 = 0: 



rlO 



8 (e Sl1 - I) 2 
5 P + P 



i 6 ; 



or 



5 Ax + lie 



5t 



hr 



h+ + h- 



-(lnu) x + 



br 



h+ + h- 



-(lntt)s 



5(Ax) 2 + 4(1 - e~ 6r ) (\nu- e 6r In a 



8 (e 6r - If 
6 h+ + h- 



[(lnw)x - (lnu) 



(47) 

As in the continuous case, there is a reduction in the difference case. When we 
consider an invariant solution, we have the reduction of the equation-grid system 
to a system of ordinary difference equations. One being the difference model for 
the considered equation, the other for the evolution of the grid. 

Let us find the solution of the difference model ()47)1 which is invariant with 
respect to the operator 



uexp 



Se St x 2 



2aX 2 + A 3 , a = const. 
Ax 



(48) 



x 



Jt 



and t are all the invariants with 



a + e St 4/ ' \e St {e Sr -l) a + e^ 
respect to (14811 . Therefore we will seek an invariant solution in the form: 



u(x, t) 
Ax 



exp 



e St( e 5r 



5e St x 2 
a + e 6t 4 

+ g(t). 



Jit)- 



X 



a + e 



St 



Substituting this form of the solution into (J47|) . we obtain a system of ordinary 
difference equations to determine fit) and git): 



f(t + T)-e Sr f(t) 
e Srt e 5r _ i) 



Jt 



• a + e' 



st- 



git) 



o. 
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The solution of this system yields the solution of the difference equation (JJJJ): 



and the grid 

e St + a 

x = x . 

1 + a 

Here x = x\ = Xi(tj) and t = tj. For t — the grid can be arbitrary, but if a 
regular grid is used the grid will be regular on every time layer. 

The obtained solution is the solution of the Cauchy problem with initial con- 
ditions: 

6e St x 2 



M (*,0) = exp /(0)-— - M - 



a + e L 

4. Invariant discrete version of the linear heat equation 

The linear heat transfer equation 

u t = u xx (49) 

admits a 6-parameter Lie symmetry group of point transformations, correspond- 
ing to the following infinitesimal operators: 

d d n d d 

Xi = -; X 2 = — ; X 3 = 2t—-xu—; 
ot ox ox Ou 

X 4 = 2t^- + x^-; X 5 = At 2 ^- + Atx^- - (x 2 + 2t)u^-; X 6 = u^-; 
ot ox ot ox ou ou 

(50) 

and an infinite-dimensional symmetry obtained from the linearity of the equation 
(Hi: 

d 

X* = a(x,t) — , 

where a(t,x) in an arbitrary solution of the equation (J49|) . 

Now we are in a position to show that the invariant difference model for the 
linear heat transfer equation cannot be constructed on an orthogonal grid. The 
model 

u — u u + — 2u + M_ 

— = — h 2 — (51) 

on the regular orthogonal mesh, which is used as an invariant model in [2], (see 
also [1], p. 363), actually does not admit operators X 3 and X 5 . 



13 



Let us check for example the symmetry that is described by X3. The operator 
X 3 generates the following transformations 



x* = x + 2ta, 

r = t, ^ 2 (52) 

This transformation destroys the orthogonality of the mesh as shown in Figure 
4. 



t t 




Figure 4: Deformation of the orthogonal mesh. 

The transformation (J52j) transforms the finite-difference equation (j51j) into 
the following: 

— 2 _T 7 

ue Ta — u u + e — 2u + u_e 

r = h* ' (53) 

which explicitly depends on a group parameter a. The first differential approxi- 
mation of the equation 

ut = u xx — Aau x + 2ua 2 + 0(r + h) 

shows explicitly the absence of invariance for the equation (JSTj) on an orthogonal 
mesh. 

Consequently, we have to construct a difference model for ()49|) on a moving 
mesh. With help of the differential invariants in the space (t, dt, dx, du): 

J l = dX+ ,lj dt i J 2 = - + ~\(-^ + 4)dt; 

dt 1 ' 1 u A dt \ u u z J 

we can represent the heat equation (}4T?|) as the system: 

J 1 = 0; 
J 2 = 0: 
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or 



dx 
~dt 

du 
~dt 



-2^; 
u 



2 U \ 
u 



(54) 



The system (|54|) inherits the set of operators Xi-Xq, X* . Implying that it entirely 
inherits the Lie symmetry group admitted by the linear heat equation (|49|) . 

For the difference modeling of the system (J54"j) we need the whole set of dif- 
ference invariants of the Lie symmetry group (joTIjl in the difference space, corre- 
sponding to the chosen stencil (t, t, x, x, h + , h~, h + , h~ ,u, u, u + , w_, u+, w_): 



h+ 



h + 

Axh + 

T 

Axh + 

T 



T 



1/2 



U 



T 



+ 



+ 



h+ + h- 

h+ 2 
h+ + h- 

2h + 
h+ + h- 

2h + 
h+ + h- 



h + u 



h + u 
1 



exp 



'1 (Ax) 2 
4 r 



1 



-ln- 



u. 



n u 
h~ u 



' h~ j u + h + m_ n 
, h + u hr u 

- — In— - — ln — 

h+ u hr u 



Approximating the system (J54|) by these invariants as it was done for the 
semilinear heat equation, we obtain a system of difference evolution equations. 
As an example, we present here an invariant difference model that has explicit 
equations for the solution and the trajectory of x: 



Ax 



It 



I ex p 



hr + hr 

1 (Ax) 2 

2 r 



h u + 

TT ln— 

n + u 



h + U- 

7— ln — 

h u 



(55) 



1 - 



At 
h+ + h~ 



1 . U + 1 , M_ 
T- ln — + TT ln — 

hr u h u 



5. Example of an exact solution. 

Let us find the solution of the difference model for the heat equation that is 
invariant with respect to the operator 

2aX 2 + X 3 , a = const. (56) 
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The operator (jo^j) has three invariants: t and the expressions uexp( i ^_ a ^ ] 
and (— — 7T-). So, we will seek the invariant solution in the form: 



t+a' 

X 2 



u(x,t) = /(t)exp 



Ax , , x 

— = 9(t) 



4(t + a) ) ' 



r t + a 

Substituting this form of the solution to the system (|55J) . we obtain ordinary 
difference equations for f(t) and g{t): 



»(() = o. 

Solving this system, we find the solution of the difference equation 

. exp ; 

t + a) V 4(t + a) 

and the solution for evolution of a grid: 



(57) 



u(x,t) = /(0) ( ) exp 



a 



The obtained solution is the solution of the Cauchy problem with the invariant 
initial condition: 



u(x, 0) = /(0) exp 



x 2 



4a y 

If a = 0, the fundamental solution of the heat equation 

i \ 1/2 / 2 



«(x,0 = C7l-l exp(--l. (58) 
is a solution of the difference model. This solution holds on the grid: 

. T 

Ax = —x. 
t 

In all cases listed above the difference mesh is arbitrary at the initial point, 
t = 0. In this case it will not be uniform on other time layers. If the grid is 
uniform in x-direction at t — (h + = h_ = h), the steps of the grid in the 
x-direction will be equal each other on every time layer, but differ from steps on 
the previous time layer. 
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It is necessary to mention that the obtained difference invariant solution is the 
solution of the corresponding differential equation that is invariant with respect 
to the operator (}56|) . As in the differential case the above reduction procedure 
can be applied for every subalgebra of the algebra (f5U|) . and then one obtain 
different moving meshes which are self-adaptive to every solution. 

Thus, the above difference models inherit both the groups of the differential 
equations and the potential to be integrated on a subgroup. 

6. Numerical calculations 

Here we do not discuss the questions of stability and convergence of the developed 
schemes. These are hard questions for nonlinear schemes but one of the ways to 
check them is by computing the numerical solutions to the exact solutions of 
the original differential equations. Below we present the results of numerical 
calculation of the invariant solution (p)H|) by means of the invariant model ()55|) . 
It is necessary to note that the calculations were not done for the equations (JoTj) 
reduced on the subgroup, but for the nonstationary equations (J57)|) . Initial data 
correspond to the solution with t = 10. In Fig. 5 we present the evolution 
of u from invariant initial data by the invariant scheme (|55j) . 




-20 -15 -10 -5 5 10 15 20 



Figure 5: Solution of the invariant model. 

Let us note that the difference model (|55|) gives us practically the exact solu- 
tion for equation (|49jl . There is only round-off error of computations. In the Fig. 
6 the evolution of the grid in the plane (t, x) for the calculation of the solution 
(J5HJ) is shown. 

The same calculations for the difference equation (joTj) on the orthogonal grid 
(as that in Fig. 1) gives a similar picture (see Fig. 7). In this case the solution 
does not coincide with exact solution of the equation fi9|) . The difference between 
exact solution and numerical solution is shown in Fig. 8. 

We should note that for the calculation on the model we defined the 
boundary values u on the moving ends of the space interval. For the difference 
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Figure 6: Evolution of the mesh. 




X 



Figure 7: Solution of the noninvariant scheme. 

equation (j51j) we defined the boundary values of u on the ends of the fixed or- 
thogonal grid in accordance with the same solution. The comparison of the two 
different models shows that for the invariant model even on the decreasing num- 
ber of the points of the grid on the initial space interval we have greater accuracy 
than for the noninvariant one. 

Remarks 
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Figure 8: The error of the noninvariant scheme. 



Following the above technique for the Burgers equation for the potential 

1 



w t + -w x 



w. 



we obtain the finite-difference model for this equation on a moving mesh 



Ax 



T- 



h w x + h + w . 

h h ' 



h+ + h 
Ax 2 



exp \ w — w 



2r 



1 + T W . 



w , 



where w x . 

h " + "' i r h 

It is interesting to note that this model is connected with model (I55j) by the 
same Hopf transformation 

w = —2 lnw 

as their continuous counterparts. 

It is important to notice that in all cases the moving in (x, t)- plane meshes 
can be stopped by the new coordinates of Lagrange's type with one additional 
dependent variable (for involving those coordinates see, for example, [14]). 
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